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In this paper, we investigate the microscopic dynamics of a polymer of length N translocating 
through a narrow pore. Characterization of its purportedly anomalous dynamics has so far re- 
mained incomplete. We show that the polymer dynamics is anomalous until the Rouse time 
tr ~ 7V 1+2y , with a mean square displacement through the pore consistent with t( 1 +")/( 1 + 2 ^) , 
with v Ri 0.588 the Flory exponent. This is shown to be directly related to a decay in time of 
the excess monomer density near the pore as t - ' 1+ ") A 1 + 21 ') exp(— t/r^). Beyond the Rouse time 
translocation becomes diffusive. In consequence of this, the dwell-time r<j, the time a translocating 
polymer typically spends within the pore, scales as 7V 2 +" , in contrast to previous claims. 



I. INTRODUCTION 

Transport of molecules across cell membranes is an es- 
sential mechanism for life processes. These molecules are 
often long and flexible, and the pores in the membranes 
are too narrow to allow them to pass through as a single 
unit. In such circumstances, the passage of a molecule 
through the pore — i.e. its translocation — proceeds 
through a random process in which polymer segments 
sequentially move through the pore. DNA, RNA and 
proteins are naturally occurring long molecules (1-5) sub- 
ject to translocation in a variety of biological processes. 
Translocation is used in gene therapy (7; 8), in deliv- 
ery of drug molecules to their activation sites (9), and as 
an efficient means of single molecule sequencing of DNA 
and RNA (6). Understandably, the process of transloca- 
tion has been an active topic of current research: both 
because it is an essential ingredient in many biological 
processes and for its relevance in practical applications. 

Translocation is a complicated process in living organ- 
isms — its dynamics may be strongly influenced by var- 
ious factors, such as the presence of chaperon molecules, 
pH value, chemical potential gradients, and assisting 
molecular motors. It has been studied empirically in 
great variety in the biological literature (10; 11). Studies 
of translocation as a biophysical process are more recent. 
In these, the polymer is simplified to a sequentially con- 
nected string of N monomers. Quantities of interest are 
the typical time scale for the polymer to leave a confin- 
ing cell or vesicle, the "escape time" (12), and the typi- 
cal time scale the polymer spends in the pore or "dwell 
time", (13) as a function of chain length N and other 
parameters like membrane thickness, membrane adsorp- 
tion, electrochemical potential gradient, etc. (14). 
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FIG. 1 Our system to study translocation in this paper. It 
consists of two cells A and B that are connected by a pore of 
diameter unity in a membrane of thickness unity. Both cells 
have the same volume V (large compared to the polymer's 
typical size). The polymer repeatedly moves back and forth 
from one cell to the other through the pore. At any time, 
exactly one monomer can be within the pore. The Kuhn 
length of the polymer and the lattice spacing are also set 
to unity. Polymers can be in three different states (i) state 
A: all monomers are in cell A; (ii) state T (threaded): some 
monomers are in cell A and the rest in cell B; (iii) state B: 
all monomers are in cell B. The dwell time is defined as 
the pore-blockade time in experiments, i.e., as how long the 
polymer spends in state T during a translocation event. 



These have been measured directly in numerous exper- 
iments (16). Experimentally, the most studied quantity 
is the dwell time t<2, i.e., the pore blockade time for a 
translocation event. For theoretical descriptions of t<j, 
during the last decade a number of mean-field type theo- 
ries (12-14) have been proposed, in which translocation is 
described by a Fokker-Planck equation for first-passage 
over an entropic barrier in terms of a single "reaction 
coordinate" s. Here s is the number of the monomer 
threaded at the pore (s = 1,...,N). These theories 
apply under the assumption that translocation is slower 
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than the equilibration time-scale of the entire polymer, 
which is likely for high pore friction. In Rcf. (17), this 
assumption was questioned, and the authors found that 
for a self- avoiding polymer performing Rouse dynamics, 
Td > tr, the Rouse time. Using simulation data in 2D, 
they suggested that the inequality may actually be an 
equality, i.e., Td ~ tr ~ N 1+2v ~ TV 2 - 18 . This suggestion 
was numerically confirmed in 2D in Ref. (18). However, 
in a publication due to two of us, Td in 3D was numer- 
ically found to scale as - TV 2 - 40±0 05 (15). Additionally, 
in a recent publication (21) r d was numerically found to 
scale as JV 2 52±0 04 in three dimensions [a discussion on 
the theory of Rcf. (21) appears at the end of Sec. IV]. 

Amid all the above results on r d mutually differing by 
~ O(N - 2 ), the only consensus that survives is that > 
tr (15; 17). Simulation results alone cannot determine 
the scaling of r^: different groups use different polymer 
models with widely different criteria for convergence for 
scaling results, and as a consequence, settling differences 
of ~ 0(A 2 ) in 0(tr), is extremely delicate. 

An alternative approach that can potentially settle the 
issue of Td scaling with N is to analyze the dynamics of 
translocation at a microscopic level. Indeed, the lower 
limit tr for Td implies that the dynamics of translocation 
is anomalous (17). We know of only two published stud- 
ies on the anomalous dynamics of translocation, both us- 
ing a fractional Fokker-Planck equation (FFPE) (20; 21). 
However, whether the assumptions underlying a FFPE 
apply for polymer translocation are not clear. Addition- 
ally, none of the studies used FFPE for the purpose of 
determining the scaling of Td- In view of the above, such 
a potential clearly has not been thoroughly exploited. 

The purpose of this paper is to report the characteris- 
tics of the anomalous dynamics of translocation, derived 
from the microscopic dynamics of the polymer, and the 
scaling of Td obtained therefrom. Translocation proceeds 
via the exchange of monomers through the pore: imagine 
a situation when a monomer from the left of the mem- 
brane translocates to the right. This process increases 
the monomer density in the right neighbourhood of the 
pore, and simultaneously reduces the monomer density 
in the left neighbourhood of the pore. The local en- 
hancement in the monomer density on the right of the 
pore takes a finite time to dissipate away from the mem- 
brane along the backbone of the polymer (similarly for 
replenishing monomer density on the left neighbourhood 
of the pore). The imbalance in the monomer densities 
between the two local neighbourhoods of the pore dur- 
ing this time implies that there is an enhanced chance 
of the translocated monomer to return to the left of the 
membrane, thereby giving rise to memory effects, and 
consequently, rendering the translocation dynamics sub- 
diffusive. More quantitatively, the excess monomer den- 
sity (or the lack of it) in the vicinity of the pore manifests 
itself in reduced (or increased) chain tension around the 
pore, creating an imbalance of chain tension across the 
pore (we note here that the chain tension at the pore 
acts as monomeric chemical potential, and from now on 



we use both terms interchangeably). We use well-known 
laws of polymer physics to show that in time the im- 
balance in the chain tension across the pore relaxes as 
t -(i+u)/(i+2u) exp(-£/Tij) (22). This results in translo- 
cation dynamics being subdiffusivc for t < tr, with the 
mean-square displacement (As 2 (t)) of the reaction co- 
ordinate s(t) increasing as £( 1+I/ )/( 1 + 2l/ ); and diffusive for 
t > tr. With ^(As 2 (T d )) ~ N, this leads to T d ~ N 2+l/ . 

This paper is divided in four sections. In Sec. II we 
detail the dynamics of our polymer model, and outline its 
implications on physical processes including equilibration 
of phantom and self-avoiding polymers. In Sec. Ill we 
elaborate on a novel way of measuring the dwell time that 
allows us to obtain better statistics for large values of N. 
In Sec. IV we describe and characterize the anomalous 
dynamics of translocation, obtain the scaling of r d with 
TV and compare our results with that of Ref. (21). In 
Sec. V we end this paper with a discussion. 




FIG. 2 Illustration of the two-dimensional version of the lat- 
tice polymer model. In the upper polymer, interior monomers 
2, 4, 6, 9, 10 and 11 can either diffuse along the polymer back- 
bone, or move sideways; monomer 7 can join either 6 or 8; 
the end monomers 1 and 12 can move to any vacant nearest- 
neighbor site. In the lower polymer, interior monomers 3, 5, 6, 
10 and 11 can either diffuse along the tube, or move sideways; 
monomer 1 can move to any vacant nearest-neighbor site, and 
monomer 12 can join its neighbor 11. All other monomers are 
not mobile. 



II. OUR POLYMER MODEL 

Over the last years, we have developed a highly ef- 
ficient simulation approach of polymer dynamics. This 
approach is made possible via a lattice polymer model, 
based on Rubinstein's repton model (30) for a single 
reptating polymer, with the addition of sideways moves 
(Rouse dynamics) and entanglement. A detailed descrip- 
tion of this lattice polymer model, its computationally 
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FIG. 3 Collapse of c(t) for different values of N, showing that the equilibration times for phantom and self-avoiding polymers 
scale as -/V 2 and N 1+2v respectively. Here N is the polymer length, (a) Data for phantom polymers, simulations were run with 
the same definition of time for all values of N; to achieve the data collapse, the times were then scaled by a factor TV 2 , (b) 
Data for self-avoiding polymers, simulations were run with the same definition of time for all values of N; to achieve the data 
collapse, the times were scaled by a factor N 1+2v (N 2 ' 2 to be precise). Note that the units of time are arbitrary and clearly 
not important for the scaling of polymer equilibration times. 



efficient implementation and a study of some of its prop- 
erties and applications can be found in Refs. (23; 24). 

In this model, polymers consist of a sequential chain 
of monomers, living on a FCC lattice. Monomers adja- 
cent in the string are located either in the same, or in 
neighboring lattice sites. The polymers are self-avoiding: 
multiple occupation of lattice sites is not allowed, ex- 
cept for a set of adjacent monomers. The polymers move 
through a sequence of random single-monomer hops to 
neighboring lattice sites. These hops can be along the 
contour of the polymer, thus explicitly providing repta- 
tion dynamics. They can also change the contour "side- 
ways" , providing Rouse dynamics. Time in this polymer 
model is measured in terms of the number of attempted 
reptation moves. A two-dimensional version of our three- 
dimensional model is illustrated in fig. 2. 



A. Influence of the accelerated reptation moves on 
polymer dynamics 

From our experience with the model we already know 
that the dynamical properties are rather insensitive to 
the ratio of the rates for Rouse vs. reptation moves (i.e., 
moves that alter the polymer contour vs. moves that 
only redistribute the stored length along the backbone). 
Since the computational efficiency of the latter kind of 
moves is at least an order of magnitude higher, we exploit 
this relative insensitivity by attempting reptation moves 
q times more often than Rouse moves; typical values are 
q = 1, 5 or 10, which correspond to a comparable amount 
of computational effort invested in both kinds of moves. 
Certainly, the interplay between the two kinds of moves 



is rather intricate (25). Recent work by Drzewinski and 
van Leeuwen on a related lattice polymer model (26) pro- 
vides evidence that the dynamics is governed by Nq~ x / 2 , 
supporting our experience that, provided the polymers 
are sizable, one can boost one mechanism over the other 
quite a bit (even up to q ~ N 2 ) before the polymer dy- 
namics changes qualitatively. In order to further check 
the trustworthiness of this model, we use it to study the 
equilibration properties of polymers with one end teth- 
ered to a fixed infinite wall (this problem relates rather 
directly to that of a translocating polymer: for a given 
monomer threaded into the pore, the two segments of the 
polymer on two sides of the membrane behave simply as 
two independent polymer chains; see Fig. 1). This par- 
ticular problem, wherein polymer chains (of length N) 
undergo pure Rouse dynamics (i.e., no additional repta- 
tion moves) is a well-studied one: the equilibration time 
is known to scale as N 1+2l/ for self-avoiding polymers 
and as TV 2 for phantom polymers. To reproduce these 
results with our model we denote the vector distance of 
the free end of the polymer w.r.t. the tethered end at 
time t by e(t), and define the correlation coefficient for 
the end-to-end vector as 

c(t) = (e(t)-e(0))-(e(f)He(0)) 

v /(e 2 (i)-(e(i)} 2 )(e 2 (0)-(e(0)) 2 )- 

The angular brackets in Eq. (1) denote averaging in 
equilibrium. The c(t) quantities appearing in Fig. 3 
have been obtained by the following procedure: we first 
obtain the time correlation coefficients c(t) for 32 inde- 
pendent polymers, and c(t) is a further arithmatic mean 
of the corresponding 32 different time correlation coeffi- 
cients. For different values of N we measure c(t) for both 
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FIG. 4 A typical translocation process of the polymer as for our system the polymers move repeatedly back and forth between 
cells A and B. See Fig. 1 and text for the definitions of the states A, B, T, M, T' and M'. 



self- avoiding and phantom polymers. When we scale the 
units of time by factors of A^ 2 and N 2 - 2 respectively for 
phantom and self-avoiding polymers, the c(t) vs. t curves 
collapse on top of each other. This is shown in Fig. 3. 
Note here that 1 + 2v = 2.175, which is sufficiently close 
to 2.2, and in simulations of self-avoiding polymers [Fig. 
2(b)] we cannot differentiate between l + 2v and 2.2. 



III. TRANSLOCATION, DWELL AND UNTHREADING 
TIMES 

Our translocation simulations are carried out only for 
self-avoiding polymers. For long polymers, full translo- 
cation simulations, i.e., having started in one of the cells 
(Fig. 1), finding the pore and distorting their shapes to 
enter the pore to finally translocate to the other side 
of the membrane are usually very slow. An uncoiled 
state, which the polymer has to get into in order to pass 
through the pore is entropically unfavourable, and there- 
fore, obtaining good translocation statistics for translo- 
cation events is a time-consuming process. To overcome 
this difficulty, in Rcf. (15), we introduced three different 
time scales associated with translocation events: translo- 
cation time, dwell time and unthreading time. For the 
rest of this section, we refer the reader to Fig. 1. 



A. Translocation and dwell times 

In states A and B (Fig. 1), the entire polymer is lo- 
cated in cell A, resp. B. States M and M' are defined as 
the states in which the middle monomer is located exactly 
halfway between both cells. Finally, states T and T' are 
the complementary to the previous states: the polymer is 
threaded, but the middle monomer is not in the middle of 
the pore. The finer distinction between states M and T, 
resp. M' and T' is that in the first case, the polymer is on 
its way from state A to B or vice versa, while in the sec- 
ond case it originates in state A or B and returns to the 
same state. The translocation process in our simulations 



can then be characterized by the sequence of these states 
in time (Fig. 4). In this formulation, the dwell time r d 
is the time that the polymer spends in states T, while 
the translocation time Tt is the time starting at the first 
instant the polymer reaches state A after leaving state 
B, until it reaches state B again. As found in Ref. (15), 
Td and r t are related to each other by the relation 

T t = VT d N 1 -*-**, (2) 

where 7 = 1.1601, 71 = 0.68 and V is the volume of cell 
A or B (see Fig. 1). 

B. Unthreading time and its relation to dwell time 

The unthreading time t u is the average time that either 
state A or B is reached from state M (not excluding pos- 
sible recurrences of state M). Notice that in the absence 
of a driving field, on average, the times to unthread to 
state A equals the the times to unthread to state B, due 
to symmetry. The advantage of introducing unthreading 
time is that when one measures the unthreading times, 
the polymer is at the state of lowest entropy at the start 
of the simulation and therefore simulations are fast and 
one can obtain good statistics on unthreading times fairly 
quickly. Additionally, the dwell and unthreading times 
are related to each other, as outlined below, and using 
this relation one is able to reach large values of N for 
obtaining the scaling of the dwell time. 

The main point to note that the dwell time can be 
decomposed into three parts as 

Td = 7"A— >M + 7"MR->M + T~m^B , (3) 

whereas mean unthreading time can be decomposed into 
two parts as 

t u = tm^m + Tm->b ■ (4) 

Here ta^m, tm^m and tm-+b respectively are the mean 
first passage time to reach state M from state A, mean 
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time between the first occurance at state M and the last 
occurance of state M with possible reoccurances of state 
M, and the mean first passage time to reach state B from 
state M without any reoccurance of state M. Since on av- 
erage ta— >m = t~m^b due to symmetry, and all quantities 
on the r.h.s. of Eqs. (3) and (4) are strictly positive, we 
arrive at the inequality 



t u < r d < 2t u . 



(5) 



Since Eq. (5) is independent of polymer lengths, on av- 
erage, the dwell time scales with N in the same way as 
the unthreading time (27). 



IV. CHARACTERIZATION OF THE ANOMALOUS 
DYNAMICS OF TRANSLOCATION AND ITS RELATION 

TO r d 

The reaction coordinate [the monomer number s(t'), 
which is occupying the pore at time t'] a convenient choice 
for the description of the microscopic movements of the 
translocating polymer, since the important time-scales 
for translocation can be obtained from the time evolu- 
tion of the reaction coordinate s(t') as shown below. To 
delve deeper into its temporal behavior, we determine 
Pn,t(si, S2,t), the probability distribution that at time t 
a polymer of length N is in a configuration for which 
monomer si is threaded into the pore, and the poly- 
mer evolves at time t' + 1 into a configuration in which 
monomer S2 is threaded into the pore. The subscript r 
denotes our parametrization for the polymer movement, 
as it determines the frequency of attempted reptation 
moves and the sideways (or Rouse) moves of the poly- 
mer. See Sec. II. A for details. 

To maintain consistency, all simulation results re- 
ported here in this section are for q ~ 10, so we drop 
q from all notations from here on. This value for q is 
used in view of our experience with the polymer code: 
q = 10 yields the fastest runtime of our code; this was 
the same value used in our earlier work (15). As dis- 
cussed in Sec. II. A, this value of q does not change any 
physics. 

First, we investigate the shape of these probability dis- 
tributions for various values of s%, S2 and t, for different 
sets of N and r values. We find that as long as the t val- 
ues are such that neither s\ nor S2 are too close to the end 
of the polymer, P/v(si, S2, i) depends only on (s2 — si), 
but the centre of the distribution is slightly shifted w.r.t. 
the starting position si by some distance that depends 
on Sx, N and q. This is illustrated in Fig. 5, where (on 
top of each other) we plot P/v(si, S2, t) for si = A/4, N/2 
and 3A/4 at t = 100 time units, for N = 400, as well as 
the Gaussian distribution. Notice in fig. 5 that the dis- 
tribution P/v(si, S2,t) differs slightly from Gaussian (the 
parameter for the Gaussian is calculated by least-square 
optimization) . We find that this difference decreases with 
increasing values of t (not shown in this paper). 



We now define the mean (\s2 — s\\){s\,t) and the vari- 
ance (As 2 (si,t)) of the distribution Pn(si, S2,t) as 

([»2 - si])(si,t) = J ds 2 PN(s l7 s 2 ,t) (s 2 - si); 



<As 2 (*i,i)> 



ds 2 PN(si,S2,t) X 



{(S2-Sif- {(S2-Si)( si ,t)f}, (6) 



where the quantities within parenthesis on the l.h.s. of 
Eq. (6) denote the functional dependencies of the mean 
and the standard deviation of (s 2 — We also note 
here that we have checked the skewness of Pjv(si, S2, t), 
which we found to be zero within our numerical abilities, 
indicating that Pn(si, S2, t) is symmetric in (s2 — Si). 
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FIG. 5 P N (si,s 2 ,t) for TV = 400, at t = 100. Note the 
data collapse when plotted as a function of (s2 — si). The 
distribution differs slightly from Gaussian. 

In principle, both the mean and the standard deviation 
of (s 2 — si) can be used to obtain the scaling of with N, 
but the advantage of using (As 2 (si,t)) for this purpose is 
that, as shown in Fig. 5, it is independent of s\, so from 
now on, we drop s\ from its argument. Since unthreading 
process starts at si = N/2, the scaling of with TV is 
easily obtained by using the relation 



<A S 2 (r d )) ~ N 2 



(7) 



in combination with the fact that the scalings of and 
t u with N behave in the same way [inequality (5)]. Note 
here that Eq. (7) uses the fact that for an unthreading 
process the polymer only has to travel a length N/2 along 
its contour in order to leave the pore. 



A. The origin of anomalous dynamics and the relaxation of 
excess monomer density near the pore during translocation 

The key step in quantitatively formulating the anoma- 
lous dynamics of translocation is the following observa- 
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tion: a translocating polymer comprises of two polymer 
segments tethered at opposite ends of the pore that are 
able to exchange monomers between them through the 
pore; so each acts as a reservoir of monomers for the 
other. The velocity of translocation v(t) = i>(t), rep- 
resenting monomer current, responds to <fr(t), the imbal- 
ance in the monomeric chemical potential across the pore 
acting as "voltage". Simultaneously, <p(t) also adjusts in 
response to v (t). In the presence of memory effects, they 
are related to each other by <f>{t) = J Q dt'fi(t — t')v(t') via 
the memory kernel fJ,(t), which can be thought of as the 
(time-dependent) 'impedance' of the system. Supposing 
a zero-current equilibrium condition at time 0, this rela- 
tion can be inverted to obtain v(t) = J Q dt'a(t — t')(f>(t'), 
where a(t) can be thought of as the 'admittance'. In 
the Laplace transform language, ft(k) — a _1 (fc), where 
k is the Laplace variable representing inverse time. Via 
the fluctuation-dissipation theorem, they are related to 
the respective autocorrelation functions as — t') — 
{cj)(t)(j){t')} v= o and a(t - t') = (w (*)«(*') >4>=o- 

10 2 i 




1Q u I 

10 1 io 2 io 3 io 4 io 5 io 6 
t 

FIG. 6 (colour online) Simulation results for the average 
chain tension component perpendicular to the membrane 
proxied by (z' 5 ^(oo) — z^\t)) following monomer injection at 
the pore corresponding to v(t) = pS(t), with p = 10. See text 
for details. Red circles: N/2 = 50, blue circles: N/2 = 100, 
green circles: N/2 = 150, solid black line: t~ {1+v)/(1+2v) with 
v = 0.588 for self-avoiding polymers. The steeper drop at 
large times correspond to the exponential decay exp(— t/rjj) 
(we have checked this, but have not shown in this letter). 

The behaviour of fi(t) may be obtained by consider- 
ing the polymer segment on one side of the membrane 
only, say the right, with a sudden introduction of p extra 
monomers at the pore, corresponding to impulse current 
v(t) = pS(t). We then ask for the time-evolution of the 
mean response (<5$( r ) (t)} , where 5& r \t) is the shift in 
chemical potential for the right segment of the polymer 
at the pore. This means that for the translocation prob- 
lem (with both right and left segments), we would have 
(j)(t) = 5& r \t) - 6$W(i), where 5& l >(t) is the shift in 



chemical potential for the left segment at the pore due 
to an opposite input current to it. 

We now argue that this mean response, and hence 
fi(t), takes the form fj,(t) ~ t~ a exp(— t/rp>). The ter- 
minal exponential decay exp(— t/rn) is expected from 
the relaxation dynamics of the entire right segment of 
the polymer with one end tethered at the pore [see Fig. 
3(b)]. To understand the physics behind the exponent 
a, we use the well-established result for the relaxation 
time t n for n self- avoiding Rouse monomers scaling as 
t n ~ n 1+2u . Based on the expression of t n , we antic- 
ipate that by time t the extra monomers will be well 
equilibrated across the inner part of the chain up to 
rif ~ t 1 /( 1 + 2,y ) monomers from the pore, but not sig- 
nificantly further. This internally equilibrated section 
of n t + p monomers extends only r(n t ) ~ n", less than 
its equilibrated value {n t + p) , because the larger scale 
conformation has yet to adjust: the corresponding com- 
pressive force from these nt + p monomers is expected 
by standard polymer scaling (29) to follow //ffcgT) ~ 
Sr(n t )/r 2 (n t ) ~ vp/[n t r{n t )] ~ t -(l+")/(i+2v\ This 
force / must be transmitted to the membrane, through 
a combination of decreased tension at the pore and in- 
creased incidence of other membrane contacts. The frac- 
tion borne by reducing chain tension at the pore leads 
us to the inequality a > (1 + v)/(l + 2v), which is sig- 
nificantly different from (but compatible with) the value 
Q'i = 2/(1 + 2v) required to obtain ~ tr . It seems 
unlikely that the adjustment at the membrane should 
be disproportionately distributed between the chain ten- 
sion at the pore and other membrane contacts, leading 
to the expectation that the inequality above is actually 
an equality. 

We have confirmed this picture by measuring the 
impedance response through simulations. In Ref. (28), 
two of us have shown that the centre-of-mass of the first 
few monomers is an excellent proxy for chain tension at 
the pore and we assume here that this further serves as 
a proxy for (5$. Based on this idea, we track (S$^(t)) 
by measuring the distance of the average centre-of-mass 
of the first 5 monomers from the membrane, (z^ 5 '(t)), 
in response to the injection of extra monomers near the 
pore at time 0. Specifically we consider the equilibrated 
right segment of the polymer, of length N/2 — 10 (with 
one end tethered at the pore), adding p = 10 extra 
monomers at the tethered end of the right segment at 
time 0, corresponding to p = 10, bringing its length 
up to N/2. Using the proxy (z^ 5 \t)) we then track 
(6<S>( r '(t)). The clear agreement between the exponent 
obtained from the simulation results with the theoreti- 
cal prediction of a = (1 + + 2v) can be seen in 
Fig. 6. We have checked that the sharp deviation of the 
data from the power law t _ ( 1 + l/ )/( 1 + 21 ') at long times is 
due to the asymptotic exponential decay as exp(— 4/tr), 
although this is not shown in the figure. 

Having thus shown that fj,(t) ~ t~ 1 + 2 " 1 - exp(— t/rp/), we 
can expect that the translocation dynamics is anomalous 
for t < tr, in the sense that the mean-square displace- 
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ment of the monomers through the pore, (As 2 (t)) ~ t 13 
for some (3 < 1 and time t < tr, whilst beyond the Rouse 
time it becomes simply diffusive. The value = a = 
follows trivially by expressing (As 2 (t)) in terms of 
(translocative) velocity correlations (v(t)v(t')), which (by 
the Fluctuation Dissipation theorem) are given in terms 
of the time dependent admittance a(t — t'). and hence 
inversely in terms of the corresponding impedance. 




t 



FIG. 7 (colour online) Double-logarithmic plot of the mean 
squared displacement of the reaction coordinate (As 2 (t)} as 
a function of time t, for N = 100 (orange), 200 (red) and 500 
(blue). The thick black line indicates the theoretically ex- 
pected slope corresponding to (As 2 (t)) ~ t { - 1+v)/(1+2v) . The 
dashed black line corresponds to (As 2 (i)) ~ ^ 2 /( 1 + 2l ') ) w hich 
would have been the slope of the {As 2 (t)) vs. t curve in a 
double- logarithmic plot, if Ta were to scale as tr ~ N 1+2v . 

Indeed, as shown in Fig. 7, a double-logarithmic plot of 
(As 2 (t)) vs. t is consistent with (As 2 (t)) ~ t ( - 1+v ^ 1+2v \ 
The behaviour of (As 2 (t)) at short times is an artifact 
of our model: at short times reptation moves dominate, 
leading to a transport mechanism for "stored lengths" 
(30) along the polymer's contour in which individual 
units of stored length cannot pass each other. As a result, 
the dynamics of s(t), governed by the movement of stored 
length units across the pore, is equivalent to a process 
known as "single-file diffusion" on a line, characterized 
by the scaling (As 2 (t)) ~ t x l 2 (not shown here). At long 
times the polymer tails will relax, leading to (As 2 (t)) ~ t 
for t > tr. The presence of two crossovers, the first one 
from {As 2 (t)) ~ t 1 ' 2 to (As 2 (t)) ~ t(i+«0/(i+2«0 an d the 
second one from (As 2 (t)) ~ t( 1 +")/( 1 + 21 ') to (As 2 (t)) ~ t 
at t k tjj, complicates the precise numerical verification 
of the exponent (1 + v)/{l + 1v), However, as shown in 
Fig. 7, there is an extended regime in time at which the 
quantity H 1+l/)/(1+2l/) (As 2 (i)} is nearly constant. 

The subdiffusive behaviour (As 2 (t)) ~ t^r^ for t < 
tji, combined with the diffusive behaviour for t > tr 
leads to the dwell ti me scaling as Td ~ N 2+v : based on 
the criterion that \J (As 2 (t^)) ~ N. The dwell time ex- 
ponent 2 + v ~ 2.59 is in acceptable agreement with 
the two numerical results on Td in 3D as mentioned in 
the introduction of this letter, and in Table I below we 
present new high-precision simulation data in support of 
t c i ~ N 2+u , in terms of the median unthreading time. 



The unthreading time t u is defined as the time for the 
polymer to leave the pore with s(t = 0) = N/2 and the 
two polymer segments equilibrated at t = 0. Both t u and 
Td scale the same way, since t u < Td < 2t u [see Eq. (5)]. 



AT 






100 


65136 


0.434 


150 


183423 


0.428 


200 


393245 


0.436 


250 


714619 


0.445 


300 


1133948 


0.440 


400 


2369379 


0.437 


500 


4160669 


0.431 



Table I: Median unthreading time over 1,024 runs for each N. 

Our results have two main significant implications: 

(i) Even in the limit of small (or negligible) pore fric- 
tion, equilibration time scale of a polymer is smaller 
than its dwell time scale. Yet, quasi-equilibrium 
condition cannot be assumed as the starting point 
to analyze the dynamics of unbiased translocation, 
as has been done in the mean-field theories. 

(ii) Since a = (1 + v)/(l + 2v) < 1, the diffusion 
in reaction-coordinate space is anomalous. This 
means that the dynamics of the translocating poly- 
mer in terms of its reaction co-ordinate cannot be 
captured by a Fokker-Planck equation in the limit 
of small (or negligible) pore friction. 

In view of our results, a Fokker-Planck type equa- 
tion [such as a fractional Fokker-Planck equation 
(FFPE)] to describe the anomalous dynamics of 
a translocating polymer would definitely need in- 
put from the physics of polymer translocation. It 
therefore remains to be checked that the assump- 
tions underlying a FFPE does not violate the basic 
physics of a translocating polymer. 

B. Comparison of our results with the theory of Ref. (21) 

We now reflect on the theory presented in Ref. (21). 

We have defined Td as the pore-blockade time in ex- 
periments; i.e., if we define a state of the polymer with 
s(t) = as '0' (polymer just detached from the pore on 
one side), and with s(t) — N as 'N', then Td is the first 
passage time required to travel from state to state N 
without possible reoccurances of state 0. In Ref. (21), the 
authors attach a bead at the s = end of the polymer, 
preventing it from leaving the pore; and their translo- 
cation time (r v hereafter) is defined as the first passage 
time required to travel from state to state N with re- 
occurances of state 0. This leads them to express t v in 
terms of the free energy barrier that the polymer en- 
counters on its way from state to s = N/2, where the 
polymer's configurational entropy is the lowest. Below 
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we settle the differences between r v of Ref. (21) and our 

Consider the case where we attach a bead at s = and 
another at s = N, preventing it from leaving the pore. 
Its dynamics is then given by the sequence of states, e.g., 



leading to Td ~ N 1+2u (whose resemblance to the Rouse 
time is a coincidence — note that with hydrodynamics 
Rouse time loses all relevance). These results, however, 
do need to be verified by simulations incorporating hy- 
drodynamics. 



.N x m x x'O x'm'x'm'x'O x'OxmxmxmxNx'N. 



References 



where the corresponding times taken (t„ and Td) are indi- 
cated. At state x and x' the polymer can have all values of 
s except 0, N/2 and N; and at states m and m', s = N/2. 
The notational distinction between primed and unprimed 
states is that a primed state can occur only between two 
consecutive states 0, or between two consecutive states 
N, while an unprimed state occurs only between state 
and state N. A probability argument then leads us to 



/* (1 + /m) 



Td Px + Pm (Pm + Pm')/m(l + /x) 



(8) 



where p m , p m i and p x are the probabilities of the corre- 
sponding states, / m = Pm/fW and f x = p m /p x . Since 
the partition sum of a polymer of length n with one end 
tethered on a membrane is given by Z n ~ A" n 7l_1 with 
A a non-universal constant and 71 = 0.68 (31), we have 

l/N. Similarly, 



Pm + Pm 

/ x ~ l/N (15). Finally, f m « 1 (32) yields t v ~ T d . 

We have thus shown that the free energy barrier does 
not play a role for t v , implying that the theoretical ex- 
pression for t v in Ref. (21) cannot be correct. The nu- 
merical result t v ~ 7V 2 ' 52±0 - 04 in Ref. (21), however, 
confirms our theoretical expression Td ~ N 2+v . 



V. CONCLUSION 

To conclude, we have shown that for the swollen Rouse 
chain, translocation is sub-diffusive up the configura- 
tional relaxation time (Rouse time) of the molecule, after 
which it has a further Fickian regime before the longer 
dwell time is exhausted: the mean square displacement 
along the chain is described by (As 2 (t)} ~ t( 1 + I ^)/( 1 + 2 ^) 
up to t ~ N 1+2v , after which (As 2 (i)} ~ t. Consequently, 
the mean dwell time scales as r<j ~ N 2+u . 

In future work, we will study the role of hydrodynam- 
ics. Rouse friction may be an appropriate model for the 
dynamics of long biopolymers in the environment within 
living cells, if it is sufficiently gel-like to support screened 
hydrodynamics on the timescale of their configurational 
relaxation. However, we should also ask what is expected 
in the other extreme of unscreened (Zimm) hydrodynam- 
ics. For our theoretical discussion the key difference is 
that, instead of the Rouse time tr, in the Zimm case the 
configurational relaxation times scale with N according 
to T Zimm ~ N^ u in 3D, which upon substitution into 
our earlier argument would gives the lower bound value 
a = (l + v)/(3v) for the time exponent of the impedance, 



[1] B. Dreiseikelmann, Microbiol. Rev. 58, 293 (1994). 
[2] J. P. Henry et al, J. Membr. Biol. 112, 139 (1989). 
[3] J. Akimaru et al, PNAS USA 88, 6545 (1991). 
[4] D. Goerlich and T. A. Rappaport, Cell 75, 615 (1993). 
[5] G. Schatz and B. Dobberstein, Science 271, 1519 (1996). 
[6] J. J. Nakane, M. Akeson, A. Marziali, J. Phys.: Cond. 

Mat. 15, R1365 (2003). 
[7] I. Szabo et al. J. Biol. Chem. 272, 25275 (1997). 
[8] B. Hanss et al, PNAS USA 95, 1921 (1998). 
[9] Yun-Long Tseng et al, Molecular Pharm. 62, 864 (2002). 
[10] W. T. Wickner and H. F. Lodisch, Science 230, 400 
(1995). 

[11] S. M. Simon and G. Blobel, Cell 65, 1 (1991); D. Goerlich 
and I. W. Mattaj, Science 271, 1513 (1996); K. Verner 
and G. Schatz, Science 241, 1307 (1988). 

[12] P. J. Sung and W. Park, Phys. Rev. E 57, 730 (1998); 
M. Muthukumar, Phys. Rev. Lett. 86, 3188 (2001). 

[13] P. J. Sung and W. Park, Phys. Rev. Lett. 77, 783 (1996); 
M. Muthukumar, J. Chem. Phys. Ill, 10371 (1999). 

[14] D. K. Lubensky and D. R. Nelson, Biophys. J. 77, 1824 
(1999); P. J. Park and W. Sung, J. Chem. Phys. 108, 
3013 (1998); E. Slonkina and A. B. Kolomeisky, J. Chem. 
Phys. 118, 7112 (2003). 

[15] J. K. Wolterink, G. T. Barkema and D. Panja, Phys. 
Rev. Lett. 96, 208301 (2006). 

[16] J. Kasianowicz et al, PNAS USA 93, 13770 (1996); E. 
Henrickson et al, Phys. Rev. Lett. 85, 3057 (2000); A 
Meller et al, Phys. Rev. Lett. 86, 3435 (2001); M. Akeson 
et al, Biophys. J. 77, 3227 (1999); A. Meller et al, PNAS 
USA 97, 1079 (2000); A. J. Storm et al, Nanoletters 5, 
1193 (2005). 

[17] J. Chuang et al, Phys. Rev. E 65, 011802 (2002); Y. 
Kantor and M. Kardar, Phys. Rev. E 69, 021806 (2004). 

[18] I. Huopaniemi, K. Luo, T. Ala-Nissila, S.-C. Ying, J. 
Chem. Phys. 125 124901 (2006). 

[19] The exponent v is also known as the Flory exponent and 
sometimes also referred to as the "swelling exponent" . 

[20] R. Metzler and J. Klafter, Biophys. J. 85 2776 (2003). 

[21] J. L. A. Dubbeldam, A. Milchev, V. G. Rostianshvili and 
T. A. Vilgis, e-print archive cond-mat/070166. 

[22] Strictly speaking, tr in this expression should be re- 
placed by the characteristic equilibration time of a teth- 
ered polymer with length of O(N); since both scale as 
we use tr here, favouring notational simplicity. 

[23] A. van Heukelum and G. T. Barkema, J. Chem. Phys. 
119, 8197 (2003). 

[24] A. van Heukelum et al, Macromol. 36, 6662 (2003); J. 
Klein Wolterink et al, Macromol. 38, 2009 (2005). 

[25] J. Klein Wolterink and G. T. Barkema, Mol. Phys. 103, 
3083 (2005). 

[26] A. Drzewinski and J.M.J, van Leeuwen, e-print archive 

cond-mat/0609281 (2006). 
[27] In Ref. (15) we overlooked the tm^m term of Eq. (3) to 

(mistakenly) deduce Td = 2t u . Nevertheless, due to the 



9 



independence of Eq. (5) on N our results for the scaling 

of Td in Ref. (15) remain unaffected. 
[28] D. Panja and G. T. Barkema, cond-mat/0706.3969. 
[29] P.-G. de Gennes, Scaling concepts in polymer physics 

(Cornell University Press, New York, 1979). 
[30] M. Rubinstein, Phys. Rev. Lett. 59, 1946 (1987); T. A. 

J. Duke, Phys. Rev. Lett. 62, 2877 (1989). 



[31] H.W. Diehla and M. Shpot, Nucl. Phys. B 528, 595 
(1998). 

[32] / m is a number slightly smaller than 1: if the polymer 
reaches s = N '/2 from state 0, due to the memory effects, 
it will have slightly higher chance to go back to state 
rather than to proceed to state N (15). 



